Stability of nodal quasi-particles in superconductors with coexisting orders 
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We establish a condition for the perturbative stability of zero energy nodal points in the quasi- 
particle spectrum of superconductors in the presence of coexisting commensurate orders. The nodes 
are found to be stable if the Hamiltonian is invariant under time reversal followed by a lattice 
translation. The principle is demonstrated with a few examples. Some experimental implications of 
various types of assumed order are discussed in the context of the cuprate superconductors. 



One of the most distinct properties of unconventional 
{e.g. d-wave) superconductors is the possible existence of 
nodal quasi-particles (QP) with a linear Dirac-like spec- 
trum in two dimensions, or lines of nodes in three dimen- 
sions. As the QP's dominate the low energy asymptotics 
of many physical properties, they are an essential feature 
of the state. A nodal QP spectrum in two dimensions has 
a linearly vanishing density of states at zero energy which 
produces several clear experimental fingerprints, e.g. a 
universal contribution to the thermal conductivityii^, a 
linear temperature dependence of the penetration depth 
at low temperatures^, and a Vi? dependence of the heat 
capacity in a magnetic field^. The spectrum can also be 
observed directly in ARPES and STM. 

The cuprate high temperature superconductors are 
known to be d-wave superconductors, and the proper- 
ties of their nodal QPs have been carefully explored in 
many experiments. Many theories have been proposed 
to account for the mechanism of superconductivity and 
the nature of the anomalous normal state. Some of these 
theories involve another type of ordering, that can ei- 
ther compete with superconductivity, coexist with it, or 
enhance it. Some examples, which we will treat ex- 
plicitly, arc antiferromagnetism, charge or spin density 
wave orders (e.g. "stripes"), and at least two forms of 
time-reversal symmetry breaking orbital antiferromag- 
netism: "d-density wave" (dDW) with an ordering vector 
q = (tt, tt)^, and "Varma loops"— which does not break 
the translational symmetry of the crystal, i.e. q = (0,0), 
but at the same time has no net orbital moment. Var- 
ious experimental and theoretical studies have provided 
evidence (sometimes conclusive, sometimes suggestive) of 
the existence of such "competing ordered states"—. 

These observations lead us to address the following 
question: how is the quasi-particle spectrum of a d-wave 
superconductor generically affected by the presence of a 
coexisting order? Since the nodal QPs are gapless, it 
is not surprising that there exist (as we shall show) cer- 
tain classes of infinitesimal perturbations that can change 
them qualitatively, either by gapping them, or by ex- 
panding the gapless locus in k-space from a point to a 
closed line (a "Fermi surface pocket" ) . 

The question of the stability of the nodes was ad- 
dressed in several previous studies, in particular in the 
presence of a general spin-orbit coupling^, in d-wave su- 
perconductor with several specific types of coexisting 



order^aiiHaSiiias and in a vortex latticei^iii^. 

In this paper, the question of the stability of the nodal 
QPs in the presence of a competing order of weak to 
moderate strength is addressed in the mean-field approx- 
imation. Specifically, we consider the quasi-particle spec- 
trum of a 2d partially filled band in the presence of a uni- 
form pairing field with d-wave symmetry and a second 
effective field, representing the competing order, which 
couples to a fermion bilinear and which we will consider 
as a "perturbation." The competing order can break any 
symmetry, such as translation symmetry or time-reversal 
symmetry; the only condition is that it is commensurate, 
i.e. the period of the ordered state is a rational multiple 
of the original lattice constant. Under these assumptions, 
we show that if the perturbation does not "nest" any pair 
of nodes (i.e., no two nodal k-points are coupled directly 
by the perturbation), and if the perturbation is invari- 
ant under time reversal or time reversal followed by a 
translation, then the nodes are stable at least until the 
strength of the perturbation exceeds a non-zero critical 
value. However, if the perturbation is not invariant un- 
der any such relative of time reversal symmetry, then the 
nodes can become gapped or can be shifted from zero en- 
ergy for any infinitesimal amount of perturbation. (Even 
where the nodes are perturbatively stable, as the pertur- 
bation increases in strength, the location of the nodes 
in k-space generally shifts. In many cases, this eventu- 
ally leads to a situation in which a pair of nodes meet or 
satisfy a more general nesting condition; then, they can 
annihilate each other leading to a gapped spectrum.) 

We now turn to the derivation of our main result. Con- 
sider a two dimensional d-wave superconductor on a lat- 
tice, whose spectrum we assume is well approximated by 
a uniform BCS mean-field Hamiltonian, Tig- Suppose 
that we add a periodic perturbation with a fundamental 
ordering wave vector G plus its harmonics. The period is 
assumed to be commensurate with the lattice, i.e. there 
exist integers N, n, and m such that NG = (^, ^). 
(It is straightforward to extend the proof to the case, 
such as a checkerboard CDW, in which there are two in- 
dependent, commensurate ordering vectors, Gi and G2, 
plus their combined harmonics.) The system is described 
by the effective Hamiltonian 



^*Uho(k)+w(k)]*k 



(1) 



For each k in the first Brillouin zone of the broken sym- 
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metry state, ho(k) and w(k) arc 2N x 2N matrices, and 
^'jj. is the 2N component spinor operator, 



k+G' ''"k+2G' ■ 



(2) 



ipl^ = (^cjj.|, c_kx j is the usual Nambu spinor, and cj^.^ cre- 
ates an electron in the Bloch state of the non-interacting 
band with wave-vector k and spin polarization a. 

Since Tio is the effective Hamiltonian of a uniform 
nodal superconductor, it is invariant under both time 
reversal, T and translation, 5r, by any lattice vector R, 
and hence it is invariant under the combined symmetry 
transformation, Tr = T5r: 



(3) 



Taking into account the fact that T is anti-unitary, and 
so satisfies = —1, it is straightforward to see that, 
with an appropriate choice of basis, 

TRi^kTR = Ar*^ (4) 

where Ar is the tridiagonal unitary matrix 



1(726 



k R 



Ar = 



i(k-|-G)-R 



(5) 



and (72 is a Pauli matrix. Moreover, by the usual argu- 
ments, it follows from the fermion anticommutation rela- 
tions, Bloch's theorem, and Kramer's theorem that /io(k) 
is a real, traceless, block diagonal matrix. Specifically, in 
the neighborhood of a nodal point, k„, 



ho 



1^ 

h;,' 



(6) 



where hQ(k)is a real, traceless, 2x2 matrix which 
asymptotically has the Dirac form, hQ(k) ~ (k— k„) • 
(vF(73 + VA(7i), and ho(k) is a real, traceless (27V — 2) x 
(2A^ — 2) matrix, {vp and va are, respectively, the Fermi 
velocity and the gap slope, and vp ■ va = 0.) 

Thus far, we have simply reproduced the usual calcu- 
lation of the quasiparticlc spectrum in an awkward basis 
with an artificially reduced first Brillouin zone obtained 
by treating G as a reciprocal lattice vector, and corre- 
spondingly we were forced to consider N times as many 
bands. The meaning of the statement that G does not 
nest the nodal points is that, for k in the neighborhood 
of k„ , all the eigenvalues of hp have a magnitude larger 
than a non-zero, positive constant. A". The key point 
is that hp is a real, traceless matrix, so its zeros have 
codimension 2, which means, generically, one can expect 
nodal points in 2D or nodal lines in 3D. 

We now proceed to demonstrate that if the perturba- 
tion, in addition to not nesting the nodal points, is in- 
variant under 7r for some R, i.e. Tj^^WTr = W, then 
the nodal point may move from k„ to a nearby point, k„. 



but it remains a nodal point. From the invariance of W 
and Eq. ([4]), it follows that 



At^w(k)AR -w(k) 



(7) 



which immediately implies that w is traceless. To make 
further progress, we express w and Ar in block form, in 
the same way as we treated hp: 



W = I // I , Ar 

u w ' 









(8) 



where w'(k) and u;"(k) are, respectively, a 2 x 2 
and a (2A^ — 2) x (2A^ — 2) Hermitian matrix, u(k) is 
a(27V - 2) X 2 matrix, and A^^ = ia2e'^-^. Since h^'(k„) 
is gapped, for a small enough perturbation, the matrix 
hg -l-w" has an inverse and hence the low energy states in 
the region of k-space near k„ can be found asymptotically 
exactly (i.e. neglecting errors that vanish in proportion 
to E/A") by diagonalizing the effective 2x2 Hamiltonian 



h;ff(k)^h^, + w'-ut (h^; + w")" 



(9) 



In particular, at a nodal point of the full Hamiltonian, 
hcff (k) must have a zero eigenvalue. It is straightforward 
to see from Eq. ([7]) that 



[A^J^hcff A^ = 0-2 hcff (72 = -hcff. 



(10) 



This, combined with the condition that it be Hermitian, 
implies that hcff is also a real, traceless matrix, hcff (k) = 
a (k) ai+b (k) (73 where a (k) and b (k) are real functions. 

The robustness of the node follows directly. The 
quasiparticlc spectrum derived from heff(k) is £'(k) = 
± ■\/a2(k) + 6^(k), which vanishes only when a (k) = 
b (k) = 0. Since we have two tuning parameters, kx 
and ky, a solution generically exists. More specifically, 
for k near k„, a(k) = va • (k — k„) -I- 6a and &(k) = 
vp ■ (k — k„) -I- 6b where 6a and 6b are small so long as 
the perturbation is weak; the node simply moves a small 
distance in k-space. The nodes must be stable for a finite 
range of strength of the perturbation, which concludes 
our proof. Conversely, for a perturbation that nests the 
nodal points, or for which no symmetry of the form of 7r 
exits, our proof breaks down, which suggests, but does 
not prove, that the nodal structure is fundamentally al- 
tered by this sort of symmetry breaking. This argument 
is related to the Wigner-von Neumann theoremi^. 

Examples - To illustrate the above principle, we explic- 
itly compute the QP spectrum of the mean-field Hamil- 
tonian of a d-wave superconductor in the presence of 
various symmetry-breaking orders that have been con- 
sidered in the context of the cuprate high temperature 
superconductors. To be explicit, we assume an underly- 
ing band-structure such that = — 2t (cos k^ + cos ky)~ 
At' cos kx cos ky — fi. The parameters used are t = 1, 
t' = —0.25 and fi = —0.9, representing a generic cuprate- 
like band structure. The gap function is taken to be 



Ak = Aq (cos fcj; — a cos ky 



(11) 
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FIG. 1: (color online.) Spectrum of a d-wave superconductor 
along the line {kx,ky) = in the presence coexisting 

SDW order {i.e. a staggered Zeeman field) of magnitude h. 
The nodal points are stable until h ~ 0.45, when two nodal 
points meet in k-space and a gap opens. 

with Aq = 0.4 and A = 1. (In some cases, we will explore 
the effect of an orthorhombic distortion, which we will 
incorporate by letting A 7^ 1.) 

As a first example, let us consider a G = (tt, tt) 
spin density wave, represented by the perturbation 
WsDW = ^ X^k cr "^"^k+G cr'^k,cr- This perturbation man- 
ifestly breaks time reversal symmetry. However, it is 
invariant under time reversal followed by a translation 
by R = ax. Diagonalizing the effective Hamiltonian nu- 
merically we find that the nodes in this case are robust. 
Upon increasing h, the nodes are shifted from their orig- 
inal position. Due to the reflection symmetry of the per- 
turbation around the tt (1, 1) direction, the nodes are con- 
strained to move along the (0, 0) to tt (1, ±1) lines. When 
h « 0.45, the nodes reach the points ±^ (1, 1), and so are 
nested by the ordering vector G. For h > 0.45, the spec- 
trum is fully gapped. The spectrum for various values of 
h is shown in Fig. [1] along the line k = ^ {k,k). 

To illustrate further the special role played by trans- 
lation symmetry, we will next consider the two com- 
binations of spin and charge density waves ("stripes"), 
shown in Fig. [5] a & b. In both cases, time-reversal 

i ^ fa ^ ^ t • f • f 
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FIG. 2: Two patterns of unidirectional spin and charge order 
("stripes") discussed in the text. The arrows represent the 
spin density, and the size of the circles represents the charge 
density. The rectangles are the unit cells. The primitive vec- 
tors are (4a, 0) and (2a, a). 




FIG. 3: (color online.) (a,b) Two patterns of spontaneous 
orbital currents on a square lattice, (a) the {n,n) d-density 
wave pattern, and (b) is a square lattice version of the (0, 0) 
order proposed in@. (c) Qualitative low energy spectrum of a 
d-wave superconductor with coexisting orbital current order 
of type (b). Note that the pair of nodes along the line kx = ky 
remain at zero energy, while the nodes along kx = —ky are 
shifted away from zero energy in opposite directions, forming 
hole-like or electron-like pockets. 



and translation symmetry are broken, such that there 
are 4 sites in the new unit cell. However, in state-a, 
time reversal followed by translation by R = 2ax, or 
by R = ay, remains an unbroken symmetry, whereas in 
state-b, no symmetry of the form of Tr, survives. (There 
do, however, remain unbroken symmetries which com- 
bine T and reflections through a plane.) We represent 
these states by a perturbation Hamiltonian of the form: 
^ = J2ra {^anT,a — haS^ ^) , where r is the Bravais lat- 
tice vector labelling a unit cell, a = 1, .., 4 is the index of 
the basis site in each unit cell, and Ur^a = J2a '^l a a^r.a,a 
and ^ = ^'Ylia'^'i.a <y^r,a.,a arc thc local charge and 
spin densities, respectively. Following the site labelling 
scheme shown in Fig. [5J we take represent the effective 
field conjugate to the spin density by the 8 component 
vector h = (/i, 0, — /i, 0). In case a, the field conjugate 
to the density is V = {0, —V,0, —V) while for case b, 
V = {0, —V,0,V). We have computed the spectrum nu- 
merically for fixed h as a function of increasing V. As 
expected on the basis of our general theorem, in case-a 
the nodal points survive until V exceeds a critical value 
of order unity. Conversely, in case-b, where no general 
theorem insures the stability of the nodes, we find that a 
gap opens for arbitrarily small V and grows as A ~ V^. 

The last example we will consider is the case of spon- 
taneous orbital current loops. Fig. [3{a,b) shows two or- 
bital current patterns: the first (a) is known as "d-density 
wave" (dDW)^ for which G = (tt, tt) , and the second (b) 
is a G = (0, 0) pattern which, from a broken symme- 
try viewpoint, is equivalent to a state defined by Varma 
on the somewhat more complex Cu-0 lattice in Q. For 
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either of these states, the perturbation Hamiltonian is 
of the form W = —iJ2rr'a Jrr'claCr'cr + H.C., where the 
eonnectivity of the eurrent network Jrr' is determined ac- 
cording to Fig. [3ja,b) , and aU the non-zero currents have 
the same magnitude J. Based on our general principle, 
we expect that in pattern (a) the superconducting nodes 
will survive at least over a finite range of J, since it is 
invariant under time reversal followed by a translation by 
R = ax. Pattern (b), on the other hand, is a G = (0, 0) 
pattern, so the nodes may be removed immediately. 

In fact, for pattern (a), Wk = at the nodal points 
(since it has d-wave symmetry), so the nodes are triv- 
ially stable. To avoid this non-generic situation, we in- 
troduce orthorhombicity of the gap function by setting 
A = 0.5 in Eq. pip . Diagonalizing the BdG Hamiltonian 
numerically for this case, we still find that the nodes ex- 
ist up to values of J of the order of the bandwidth t. 
Since pattern (b) does not break translation symmetry, 
the corresponding effective Hamiltonian can be easily di- 
agonalized. The eigenenergies are E^^^^ = ±-\/Ck + + 
2 J {sin {kxtt) — sin {kyu) + sin [{ky — k^) a]}. We see that 
even though the degeneracy at the nodes is not lifted, 
some of the nodes (the ones that lie on the line kx — --ky) 
are shifted away from zero energy: if J > 0, then the 
node a.t kj. > [k^ < 0) is shifted to positive (negative) 
energy, respectively. The low-energy spectrum is shown 
qualitatively in Fig. [Sljc). In the case of an orthorhombic 
gap function (A 7^ 1), the nodes along kx = ky are also 
shifted to finite energies by an amount proportional to 
A — 1. Whenever a node is shifted to a finite energy, a 
hole-like or electron-like pocket is formed, giving rise to 
a finite density of states at zero energy. 

Discussion - The condition established here for the per- 
turbative stability of the nodal points is a sufficient con- 
dition, not a necessary one. However, from the examples 
considered above, we see that in several cases where the 
condition is not met, the nodes are gapped or shifted 
to finite energies upon introducing an infinitesimal per- 
turbation. Therefore, we expect that in generic cases, 
this condition is more or less necessary for the existence 
of nodal points. Thus, we can classify various weakly 
ordered states according to whether they leave the QP 



nodes intact {e.g. dDW or certain types of stripes), gaps 
them {e.g. density wave order with ordering vector G 
which nests the nodal points), or moves them away from 
the Fermi energy {e.g. Varma loops). Since the exis- 
tence of nodal points has several well-defined experimen- 
tal consequences, this information can provide a non- 
trivial consistency check on the assumed occurrence of 
various forms of coexisting order in the superconducting 
phase. In particular, in situations where the experiments 
are consistent with nodal QPs, one can rule out (or at 
least give an upper bound on) competing orders which are 
not invariant under time reversal followed by translation. 
For instance, measurements of the linear in T decay of 
the superfluid density in ultra-pure crystals of the Ortho- 
II phase of YBCO^ give evidence that the nodal points 
are gapless and tied to the Fermi energy to within an ac- 
curacy of approximately IK, or in other words a fraction 
of a percent of Aq. This can probably be converted into 
a rather stringent bound on the strength of Varma loop 
order at low temperatures in this material. 

The extension of these results to the case of an incom- 
mensurate perturbation is an interesting open problem. 

A "striped superconducting" stated!, which can be 
thought of as a spontaneously developed FFLO state, was 
considered as a candidate state to explain some extremely 
anomalous transport data22. that was recently obtained 
on the stripe-ordered material La2_2;Ba2;Cu04. In the 
proposed state, the superconducting order parameter is 
modulated with wavevector G = (7r/4, 0). A calculation 
of the spectrum in this state reveals that the nodes ex- 
pand to a series of pockets along the bare Fermi surface, 
and the single particle spectral function has low energy 
weight along finite Fermi arcs. The reason that the nodes 
are not protected, despite the time reversal symmetry of 
this state, is that it cannot be viewed as a perturbed 
version of a uniform d-wave state. 
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